NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



K80-23239 


(NASA-CB-163CC5) PHOBOS AND CEIMGS: 

ANALYSIS OP SDK FACE FEATURES, EJECTA 
DYNAMICS AND A VOLATILE LOSS MECHANISM 
Pinal Report (Planetar)' Science Inst., 

Tucson, Ariz.) 26 p HC A03/MF A01 CSCL 031 G3/91 


PHOBOS AND DEIMOS : 

ANALYSIS OP SURFACE FEATURES, EJECTA 
DYNAMICS AND A VOLATILE LOSS MECHANISM 


NASW-3272 


FINAL REPORT 


MARCH 1980 


Principal Investigator: Dr. Donald R. Davis 

Planetary Science Institute 
2030 East Speedway 
Tucson, Arizona 85719 


Co-Investigators: Dr. Clark R. Chapman 

Dr. Richard Greenberg 

Dr. Stuart J. Weidenschilling 

Planetary Science Institute 


Scientific Collaborators: Dr. Peter Thomas 

Dr. Joseph Veverka 


Unclas 

15463 


Cornell University 
Ithaca, New York 


TABLE OF CONTENTS 


Title Page 

TASK I. INTERPRETATION AND ANALYSIS OF SATELLITE 
CRATER POPULATIONS 

A. Interpretation of Crater Populations 
on Phobos 

B. Phobos: Absolute Ages of Cratered 

Surfaces 

TASK II. DYNAMICAL STUDIES OF EJECTA DISTRIBUTION 
AND SATELLITE SURFACES 

A. Trajectory Dynamics and Ejecta Distri- 
bution 

B. Satellite Lifetimes 

C. Groove Formation on Phobos 

TASK III. STUDIES OF VOLATILE RELEASE FROM LOW- 
VELOCITY IMPACTS 


*. * 




IS 


Page 

1 

3-6 

7-9 

10-18 

19-22 

23-23a 

24 


TASK I . INTERPRETATION AND ANALYSIS OP SATELLITE CRATER POPULATIONS 


i | 

A. INTERPRETATION OF CRATER POPULATIONS ON PHOBOS ! 

j 

A significant question about the? crater population on Phobos is 
whether it represents a production population or an equilibrium 

population. If it were an equilibrium population, a plausible process [j 

I * 

of crater degradation and destruction would be the effects of what I j 

call "saturation equilibrium." This concept need not imply a geometric 
saturation of craters because the crater ejecta and other features of 
the cratering process serve to degrade and obliterate pre-existing craters 
before geometric saturation could be achieved. It is presently a matter 
of controversy about what fraction of geometric saturation is 
appropriate for saturation equilibrium. A variety of studies by 
Marcus, Schultz et al . , and others have demonstrated that this 

T ' ' 

fraction varies from planet to planet (indeed from unit to unit) 

i- 1 

% \ 

depending on the various physical parameters defining the cratering f j 

process (e.g., material properties, size of largest crater of the j 

planet, population index of projectiles, etc.). ij 

Several useful comparisons of the Phobos cratering data (as . j 

shown by Thomas et al. , 1979) with other crater density results j j 

■ i 

can be made. (1) Crater densities on Phobos tend toward the -2 ’! 

slope (actually closer to -1.8) on standard cumulative frequency 

plots, consistent with simple saturation equilibrium conditions. | 

H- 

(2) The densities of craters on Phobos lie far below geometric \\ 

saturation; moreover, they lie near but consistently below an ;j 

extrapolation to small diameters of Hartmann's (1973) curve 

for large lunar uplands craters. This curve has been widely re- i 

ferred to as a reference line against which to measure crater 
densities, although it is generally understood that the processes 
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governing saturation equilibrium densities of large lunar 
craters are different from those governing populations of 
different-sized craters on‘ different bodies. (3) Although 
quantitative measurements have not been reported yet, it appears 
that all degradation states of craters are represented on Phobos 
a requirement of equilibrium conditions. (4) There appear to be 
statistically significant differences in crater densities on 
different parts of Phobos (although all are with a factor of 
2 of the average density) yet all have distribution slopes close 
to the -1.8 characteristic of the composite average for all of 
Phobos. (5) The steep slope predicted by Neukum and Wise (1976) 
for primaries at small crater diameters is not recorded on 
Phobos. If the observed craters represent a production rather 
than equilibrium population, then the Neukum and Wise "standard 
curve" is not applicable to Phobos; otherwise, the population is 
in equilibrium. 

These observations are all consistent with, though not nec- 
essarily diagnostic of, saturation equilibrium. As described 
by Chapman and Jones (1977) , the equilibrium density at a given 
diameter is the ratio of the cratering rate to obliteration rate 
times the process-dependent "amount" of obliteration necessary 

to remove a crater. Although the crater density on Phobos is 

« 

variable and rather low compared to the Hartmann reference line, 
the densities are similar to, or even greater than, some small- 
crater populations on the lunar surface. Relative to the large- 
crater processes responsible for Hartmann's saturation reference 
• ' line, the low densities of small craters on the moon and on 
Phobos could be ascribed plausibly to the greater scaled dis- 


m 
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tances of ejecta redistribution for smaller craters combined 
with a possibly steeper slope for the production population at 
small sizes. 

More diagnostic of equilibrium is a diameter-independent 
mix of crater degradation states, especially the case of a pre- 
ponderance of degraded craters. This study has not quantified 
crater morphologies yet. But the qualitative appearance of high- 
resolution images of Phobos show an abundance of fresh craters 
in the sub-2 00m diameter range that appears inconsistent with 
some other examples of saturation equilibrium populations. 

(The apparent morphologies are so striking as not to be due to 
artifacts of image processing. But consideration of processing 
artifacts must be included in future quantitative studies of 
crater morphologies . ) Larger craters do exhibit a more repre- 
sentative range of morphologies. Possibly episodic blanketing 
due to formation of the larger craters yields oscillation about 
a quasi-equilibrium state in the degradation states of smaller 
craters; the relatively fresh population might mainly reflect 
recratering since the last major blanketing episode. In this 
manner, the craters on Phobos can be deemed to be in quasi- 
equilibrium with the inherently episodic process of crater 
saturation. 

m 

Thus observations (1) through (5) can be reconciled with 
a kind of saturation equilibrium population, produced by a pro- 
duction population with a somewhat steeper slope than observed. 

On the other hand, one cannot rule out the possibility that the 
. craters on Phobos represent a production function of relatively 


r.'. 
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modest slope. There is no empirical knowledge about the size 
distributions of projectiles near Mars capable of cratering 
Phobos, so there is no constraint on the slope or the production 
function other than some rather uncertain theoretical expectations 
of steeper slopes. The relatively large fraction of reasonably 
fresh craters on Phobos combined with the relatively low density 
both are consistent with a distribution not substantially modi- 
fied from production. 

If the production population slope is near -2, then the 
observed density of craters approaches the maximum density 
possible before the probability becomes high that Phobos 
would be catastrophically disrupted by a large impact. The 
details of the relevant calculations, based on the model of 
Housen et al. (1979), are given by Thomas et al. (1979) . There- 
fore a consistent interpretation for the craters on Phobos could 
be: (a) a production population of slope ^ -2; and (b) the 
surface of Phobos was renewed by a catastrophic collision sub- 
sequent to the early heavy bombardment. 

REFERENCES 
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B * PHOBOS: ABSOLUTE AGES OF CRATERED SURFACES 

Let us consider the question of the age of units on Phobos 
by comparison with adopted cratering ages for Mars. 

Hartmann et al. (1979) adopted crater production rates on Mars 
of 1 to 4 times those on the moon for the last 3.4 b.y. Impact 
velocities on Phobos are comparable to those on Mars since the 
change in potential from Phobos' s orbit to Mars adds little to 
the impacting object's velocity. The population and flux of 
impacting bodies has been the same as that for Mars so long as 
Photos has been in orbit about Mars. 

A major difficulty is scaling crater sizes between Phobos 
and Mars. The diameter of craters on two bodies scales as the 
ratio of surface gravities to the -k power. Published values 

of k range from 0.12 to 0.25. Gault and Wedekind (1977) 

0 

reported experimental values of about 0.16. If the number of 
impacting bodies varies roughly as D , the cratering rate 
correction, Z, is given by the ratio of surface gravities to the 
-2k power. How does one extrapolate k to very low values of cr? 
Reference to Housen et al . (.1979) shows that craters counted in 
this paper are dominated by gravity scaling if the surface layers 
of Phobos are very weak, but by energy scaling if Phobos is 
strong. Lunar and Martian craters are all predominantly in 
the gravity-scaling regime. 

Particularly uncertain are the relative crater sizes due 
to different strengths of target surfaces on Phobos compared with 
the moon and Mars. Using the Housen et al. parameters for "strong 
and "weak", and applying the appropriate scaling laws, we find 


8 . 


that a given projectile forms a crater 4.9 times as large on 

a weak Phobos compared with a strong Mars (using k = 0.16) and 

2.1 times as large on a strong Phobos. Since the correction 

2 

factor Z goes as D , the target strength uncertainty alone intro- 
duces more than a factor of 5 uncertainty in age for the sur- 
face of Phobos. 

If the crater production rate on Mars is twice the lunar 
cratering rate (factor of 2 uncertainty each way) then the 

Q 

best age for the composite surface of Phobos is about 3 x 10 

years. In one extreme case (weak Phobos, k = 0.25, Martian 

cratering rate- 4 times lunar) the observed crater densities may 

indicate exposure to cratering for only l/10th the nominal dura- 

8 

tion, implying an age of 3 x 10 years. The opposite extreme 
(strong Phobos, k = 0.12, Martian and Lunar cratering rates 
equal) yields 10 times the nominal exposure, implying that the 
surface of Phobos dates well back into the presumed epoch of 
early intense bombardment. 

In these comparisons, we have assumed that the lunar cratering 
production function is given by a power-law with slope -2 rather 
than by a curve having the shape observed on the lunar maria in 
these size ranges, which approaches a slope of -3.5 at some 
diameters. Lunar cratering specialists disputq whether the 
steep portion represents primary cratering or secondary cratering. 
If it is primary cratering, as argued by Neukum, then the age 
comparisons are confused. If the impacting populations on the 
different bodies have different size distributions, then the cal- 
culation is rendered meaningless. Alternatively, as discussed 






s. 


above, the Phobos/Deimos data represent equilibrium 
populations, in which case the deduced ages are lower limits, 
In summary, the surface of Phobos (and also of Deimos) 
is quite old, perhaps 3 b.y. Ages younger than 1 b.y. as 
well as older ages approaching the age of the solar system 
are also quite possible. 
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TASK XI. DYNAMICAL STUDIES OF EJECTA DISTRIBUTION AND 
SATELLITE SURFACES 

A. TRAJECTORY DYNAMICS AND EJECTA DISTRIBUTION 

The surfaces of Phobos and Deimos are both heavily cratered, 
but differ significantly in other ways. Deimos apparently 
has been blanketed with material which has partially covered 
many small or medium-sized craters and given the surface a 
rather smooth appearance. The surface of Phobos exhibits little 
evidence of blanketing and has a more battered appearance. If 
the blanketing material on Deimos is crating ejecta, it is 
surprising that the smaller satellite has been more effective 
at retaining regolith. However, the Martian satellites are 
small bodies located deep within a gravity well. As such, they 
provide an unusual dynamical environment for the ballistic 
transport of ejecta. 

To investigate dynamical effects on ejecta transport 
and retention, we have developed a computer program to study 
the dynamics of material ejected from the Martian satellites. 

Each satellite is modeled as a rotating, homogeneous, triaxial 
ellipsoid. The gravity field of Mars is included in the 
calculation. Trajectories of particles ejected from a given 
point on the surface are determined by numerically integrating 
the equations of motion. 






n. 


Prior to the current funding period# preliminary calcu- 

* 

lations were made to determine escape velocities for the 
satellites and to study the distribution of deposits pro- 
duced by non-escaping ejecta (Housen et al . , 1979, Appendix 
B.2; Housen and Davis, 1978, Appendix B k :0. Escape velocities 
were found to depend on location on the satellite surface 
and on launch direction. Also, escape velocities are much 
more variable over the surface of Phobos than on Deimos. 

The modeled ejecta azimuthal distributions on Phobos are much 
more distorted than those on Deimos due to the proximity of 
Phobos to Mars. On satellites, material is everywhere 

gravitationally bound to the surface. 

<* 

Since the start of the present funding period, we 
have used our numerical program to investigate the hypothesis 
that at least some of the extensive sets of linear features 
discovered on the surface of Phobos could be the result 
of secondary cratering from the Stickney impact as suggested 
by Head and Cintala (1979) . The appearance of the features 
. ranges from nearly continuous "grooves" to strings of discrete 
circular pits. Studies of the surface morphology of Phobos 
(Thomas, 1978) demonstrated a relationship between the large 
crater, Stickney, and the grooves. Thomas suggests that the 
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grooves are most likely manifestations of internal fracturing 
of Phobos due to the formation of Stickney. While the fracture 
theory for the origin of the grooves seems plausible, we 
should not rule out other possible modes of formation until 
detailed analysis can be made. 

In order to simulate the formation of crater chains, 
we determine the impact sites of particles ejected from 
Stickney with different velocities. Hypervelocity impact 
studies by Gault show that during a cratering event the 
ejection velocity of particles varies while the launch 
elevation angle remains nearly constant. Therefore, in the 
program, we eject particles at a given azimuth and elevation 
angle with a range of velocities. Each trajectory is inte- 
grated until it impacts the surface or escapes from the 
satellite. The program iterates the launch velocity until 
the escape speed is determined to within a few percent and 
plots the locus of impact points for each launch azimuth. 
Different launch angles are used to map the total impact 
distribution. 

Some results are shown in Figure la for Phobos at its 
present distance and spin period. The escape speeds 
for each azimuth angle are also shown. For comparison, a 
map of the observed linear features, constructed by Thomas, 

C1978), is shown in Figure lb. For the most part, the program 
produces crater chains whose impact loci only partially 
resemble those observed. One striking difference between 
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Figures la and lb is the paucity of features produced by the 
program in the region of 270°-’0 o longitude. 

However, formation of Stickney probably did not occur in 
recent times. In fact, age dating studies via crater counting 
by Thomas (1978) , indicate that Stickney was likely formed 
^3.0 Gy ago, however, there is considerably uncertainty in 
this estimate (see Appendix A. 4) This, coupled with the 
fact that Phobos is tidally evolving inward towards Mars, 
means that the grooves were formed when Phobos had a larger 
orbital semi-major axis than the current value. The range 
of observed secular accelerations coupled with the range of 
ages for Stickney means that, Phobos was likely located between 
the synchronous orbit distance and halfway between the present 
distance and the synchronous orbit distance. Calculations 
for these two limiting cases are shown in Figure 2 . Once 
again, the modeled pattern of crater chains differs from the 
observed pattern. For example, even though features are now 
produced in the 0°-270° longitude region, these chains cross 
one another in a manner not characteristic of the features on 
Phobos. 

Another important point to consider is the rotation rate 

4 

of Phobos. Estimates of the energy involved in the formation 
of Stickney suggest that Phobos was probably knocked out of 
synchronous rotation (.see Appendix A. 3) . We illustrate 
the effect of changing the rotation rate by a factor of 2 














(both spinning up and slowing down Phobos) ‘(Figure 3) . The 
model crater chains still differ in important ways from the 
observed features. 

In short, it would appear difficult to produce by 

secondary cratering, the observed pattern of grooves on 

Phobos when the satellite is located inside the current 

orbit of Deimos. However, if we consider Phobos either as a 

distant satellite of Mars, or a slowly rotating asteroid, 

a better match to the observed groove pattern is obtained 

(Figure 4) . The near-equatorial region near longitudes 90°- 

100° is depleted in chains because most material escapes 

Phobos and only a small fraction impacts in this region. 

There are two subtle differences between Figure 4. and Figure 

o 

lb. (1) the void region in Figure lb is centered on (270 - 
0°) , while the analogous region in Figure 4 is shifted east- 
ward roughly 10°-20°. (2) Chains on Phobos are seen to cross 

one another in the vicinity of 320° longitude. This overlap 
is not seen in Figure 3. These differences may be reconciled 
by spinning up Phobos somewhat. 

REFERENCES 
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Head, J.W., and M.J. Cintala (1979). PGPI Abstract, NASA 
TM 80339, 19. 
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B. SATELLITE LIFETIMES 


An extreme hypothesis to explain the differences in the 
satellite surfaces might be that Deimos was catastrophically 
disrupted by a large impact but subsequently reaccreted since 
most of its mass must have remained in orbit about Mars. 

Sweepup of a significant fraction of Deimos mass at low 
relative velocities could explain much of the blanketing of 
Deimos surface. This suggestion raises the question of 
whether it is probable that the projectile flux in the 
vicinity of Mars would have catastrophically disrupted either 
satellite since Mars formed. This calculation cannot be 
made with any certainty for we don't know whether phobos and 
Deimos are made of strong, consolidated material or if they 
are materially weak and held together mainly by gravity. 

However, an upper bound can be set on the probability of 
catastrophic disruption by assuming the satellites are 
gravitationally bound. 

The gravitational binding energy per unit mass on Phobos 
5 

Is 4.7 x 10 erg/gm which corresponds to an effective impact 
strength (the kinetic energy per unit volume necessary to 
disrupt a body) of 1.0 x 10 erg/cm , which is about a factor 
of 30 less than solid rock but 'vlOO times that for loosely bonded 
regolith. In other words, Phobos would have to be made of 
rather strong, mostly consolidated material in order for its 
material strength to dominate its gravitational binding. 
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The equivalent impact strength for a gravitationally bound 

5 3 

Deimos is 3 x 10 erg/cm , about a factor of three less than 
for Phobos. 

The ease with which the satellites are disrupted strongly 
depends on the fraction of the collisional energy that is 
converted into the kinetic energy of fragments and ejecta. 

We will consider two extreme cases, one assuming all the input 
energy is converted, while the other assumes that only 10% 
goes into ejecta energy. 

The impacting flux on the satellites is essentially the 
same as that on Mars. We estimate the probability of a 
catastrophic impact on the satellites in two steps. We calcu- 
late the crater size on Mars produced by a body carrying just 

* 

enough energy to disrupt a satellite. Then from the density 
of craters equal to or exceeding the threshold size, we calculate, 
the probability of a body hitting targets with the cross- 
section of the satellites. Table 1 summarizes the steps in this 
calculation for two limiting cases on both Phobos and Deimos. 

The collisional kinetic energy required to disrupt the 
satellites is given in Column 2 , while Column 3 lists the crater 
diameter that this impact energy would produce on Mars. The 
crater scaling algorithm of Housen et al. (.1979') 'is used with 
gravity scaling numerical coefficients, which is appropriate 

since S/pgD << 1 assuming reasonable geologic materials. Both 

7 3 5 3 

strong (S'v 3 x 10 erg/cm ) and weak (S < 1 x 10 erg/cm ) 
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materials are considered for Mars. The fourth column gives the 
cumulative number of craters having diameters equal to or 
exceeding the corresponding threshold crater sizes from 
Column 3. These data were taken from counts of Neukum and Wise 
on the intensely cratered uplands of Mars, which are thought 
to be the oldest existing geologic units on Mars. The crater 
retention age of these areas is estimated to be 3. 5-4.0 Gy f and 
could be as young as 1.8 Gy at the small crater diameter. 

If the crater populations are in equilibrium rather than pro- 
duction, then our calculation will yield upper bounds on the 
ages of the satellites. 

The expected number of collisions capable of disrupting 
gravitationally bound satellites is of order 1 or greater if all 
the collisional energy is partitioned into ejecta energy and 
about a factor of 2 less if only 10% goes into ejecta motion. 
Hence, gravitational binding alone is sufficient to account 
for the survival of the satellites against the Mars impacting 
flux over the last 4 Gy. That one satellite (Phobos) probably 
had an impact nearly large enough to disrupt it appears con- 
sistent with the above estimates of collision frequencies. 

Note that the lifetime of Deimos is not appreciably less than 
that of Phobos even though an order of magnitude less energy is 
required to disrupt it; the impacting population and smaller 
cross-section of Deimos compensate for the lower required 
energy . 

REFERENCE: Housen, K.R., L. WHkening, C.R. Chapman, and R. 

Greenberg 0.9791. ' ICarUs , in press. 


C. GROOVE FORMATION ON PIIOBOS 


The following paper by Stuart J. Weidenschilling, published in 
NATURE, describes one phase of research carried out under Task II. 
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A possible origin 

for the grooves of Phobos 

Stmrt I, WeidenschillitiR 

Planetary Science Institute. 2030 East Speedway, Suite 201 Tucson, 
Arizona 85719 


The surface of Phobos, the inner satellite of Mure, is marked by a 
series of linear depressions or grooves, These have been inter* 
preted as surface manifestations of internal fractures 1 '*. The 
origin of these fractures has been ascribed to tidal stresses 
induced by the decay of Phobos* orbit*, or to a nearly catastro* 
phic cratering event , I suggest here that the age and distribution 
of the grooves may be explained by a hybrid origin due to 
varying tidal stresses induced after a large impact altered the 
satellite’s rotation rate. 

The orbital decay hypothesis was suggested by Soter and 
Harris 5 , who noted that Phobos is subjected to increasing tidal 
stresses as it approaches Mars. As the rate of orbital decay 
increases with time, this mechanism would require the grooves 
to be extremely young features. Later imaging at higher resolu- 
tion revealed an appreciable number of impact craters 
superimposed on the grooves. Their implied age is > 10 9 yr (ref, 
2), effectively rufing out this explanation. 

Thomas et al , 5 noted that the grooves are most prominent 
near the largest crater, Stickney, and concluded that the frac- 
tures were produced by that impact. However, the relationship 
of the grooves with respect to the crater is not simple. They arc 
neither concentric nor radial to the crater; rather, they seem to 
define three sets of planes: one parallel to the equatorial plane, 
another perpendicular to the longest axis of the satellite, and a 
third intermediate between those two 1,2 , While crater counts 
show that all of the grooves have roughly the same absolute age, 
some of them show cross-cuttting which implies that they did not 
form simultaneously 1 , If the material excavated from Stickney 
went into orbit about Mars, it should have been re-accreted by 
Phobos 4 *, if uniformly distributed, it would blanket the surface to 
a depth of about 20 m. Many of the grooves exhibit relief of 
< 10 m without apparent mantling by debris 1 , suggesting that 
they were formed or renewed after the bulk of Stickney ejecta 
had settled, . 

Apparently then the existence and location of the grooves are 
related to the Stickney impact, but they formed shortly after that 
event, with orientations determined by the shape of Phobos, 
That shape is far from spherical, resembling a triaxiai ellipsoid, 
The present rotation period of Phobos is equal to its orbital 
period. The longest axis of the satellite points towards Mars; this 
orientation minimises the internal stresses due to the combina- 
tion of topography and tidal forces. These conditions almost 
certainly prevailed before the formation of Stickney, as the time 
scale for tidal despinning is <10 6 yr (ref, 5), If the 'n?j*ict 
destroyed the condition of synchronous rotation, or simply 
induced a large iibration amplitude, the changing orientation of 
Phobos in the gravitational field of Mars would subject its 
interior to periodically varying stresses, These could cause 
faulting along planes of maximum shear stress, producing the 
systems of grooves. 

A quantitative estimate shows that this mechanism is plausi- 
ble, In the stable synchronous state, the rotational energy is 
Cn 2 / 2, where C is the maximum moment of inertia and n is the 
mean motion. The criterion for destroying synchroncity, 
producing circulation rather than Iibration, is a change in rota- 
tional energy of magnitude 3(fl-A)/; 2 /2, where A and B are 
the moments of inertia about the other principal axes 6 , Model- 
ling Phobos as an ellipsoid of mass 10 19 g and semi-major axes of 
13,5, 10.8, and 9.4 km, C=*6x 10 3,1 g cm 2 , and (J3-A)/C=* 
0.25. The change in rotational energy required to break syn- 
chroneity is, coincidentally, of the same order as the total 
rotational energy in the synchronous state. At present, these 
quantities are of the order 10 23 erg. They vary as r~ } , where r is 


the orbital radius; as the orbit of Phobos is decaying, they may 
have been smaller in the past by one order of magnitude. 

Thomas estimates the energy of the Stickney event at about 
10 16 erg.This figure should be regarded as having an uncertainty 
of at least an order of magnitude. In any case, only a small 
fraction of the impact energy would need to be partitioned into 
rotational energy of Phobos to remove it from synchroneity. For 
on impact energy of 10 JI * erg, at a velocity of 11) km s 1 and an 
Impact parameter of 10 km, the angular momentum imparted is 
- l() !6 g cm’s 1 . The present rotational angular momentum of 
Phobos is -lO 2 ’ gem’s *, but varied in the past as r 1 ’.Within 
the uncertainties of the estimate, the change in rotational 
angular momentum could have been comparable to total pre- 
impact value. In that case, the spin axis could have been appre- 
ciably displaced from the axis of maximum moment of inertia, 
The resulting nutation might account for those grooves which 
are not aligned with the principal nxes. It is unlikely that the 
lesser impacts which Phobos experienced could have removed it 
from synchronous rotation; in this sense, the Stickney event was 
unique. 

The excess rotational energy dissipated within Phobos during 
despinning after the impact was much less than the impact 
heating, and would have been negligible if uniformly dis- 
tributed, However, it was presumably released along planes of 
faulting, and might have aided the escape of volatiles from the 
grooves 1 . The time scale for despinning is - l()M0 n yr (ref. 5), 
while re-accretion of ejecta orbiting Mars would take ~ 10* yr 
sref, 4), Hence, it is reasonable that the grooves show no signs of 
burial by ejecta from Stickney. 

Since Phobos shows no recent cracking in its synchronous 
state, the present stresses must be less than those produced 
during non-synehronous rotation. For a tidaily stressed elastic 
sphere 7 , the maximum tensile and shear stresses would be a few 
times 10 4 dyn cm' 2 near the surface, with the centre in 
compression at ~ tO 6 dyn cm' 2 . These values are probably over- 
estimates, as when the long axis is towards Mars, Phobos is 
nearer to hydrostatic equilibrium than a spherical body would be 
(its present shape would be almost exactly in equilibrium at 1 ,25 
times Its present distance from Mars”), Even though Phobos is 
within the classical Roche limit for a liquid satellite (as its density 
is less than that of Mars), the stresses within it arc quite small, 
However, in non-synehronous rotation, the topographic relief of 
=« 3 km would produce shear stresses of the order of a few times 
10 5 dyn cm 2 (ref. 9). The tidal stresses would be added to this 
topographic stress when the long axis of the satellite was pointed 
away from Mars, 

The region opposite Stickney shows few grooves; therefore a 
shear strength > 10* dyn cm" 5 is implied for the undisturbed 
material. This value is consistent with a carbonaceous chondrite 
composition, but a stronger material cannot be ruled out, The 
region near the crater is presumably weaker d ue to fracturing by 
the impact, and may not constrain the properties of the pristine 
material, The low strength implied is for repeated loading; the 
yield strength under unidirectional stress may be greater. More 
realistic calculations of the stress field in a tidaily distorted, 
ellipsoidal elastic body are needed to refine these strength 
estimates, and to determine whether the groove orientations are 
indeed consistent with this origin, 

If ' j grooves are simply due to impact fractures, then such 
features might be common on asteroids, as predicted by Thomas 
ej al? If a strong tidal gravity field is also necessary, as suggested 
here, then only close satellites could possess regular systems of 
grooves, The absence of grooves on Deimos could then be due to 
its greater distance from Mars and smaller size, which resulted in 
lower stresses, or lack of an impact large enough to remove it 
from synchronous rotation, The best additional candidates for 
grooves would be Amalthea, the inner satellites of Saturn, and 
possibly the larger particles in planetary rings. 

I thank C. R. Chapman, D. R. Davis and R, Greenberg for 
helpful discussions, This work was supported by NASA contract 
NASW-3272. 
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TASK III. STUDIES OF VOLATILE RELEASE FROM LOW-VELOCITY IMPACTS 


Work was initiated in reviewing data films obtained at Ames 
Research Center on impact experiments into regolith containing 
trapped gases. Analysis is in progress on the ability of sub- 
surface volatiles to displace and transport regolith, the mor- 
phology of the resulting surface features and the subsequent 
evolution of the displaced material. Applications of these 
results to the Martian Satellites and possibly to comets will be 
carried out in the next portion of this program. 


